Non-equilibrium Transport in the Anderson model of a biased Quantum Dot: 

Scattering Bethe Ansatz Phenomenology 
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We derive the transport properties of a quantum dot subject to a source-drain bias voltage at 
zero temperature and magnetic field. Using the Scattering Bethe Anstaz, a generalization of the 
traditional Thermodynamic Bethe Ansatz to open systems out of equilibrium, we derive results 
for the quantum dot occupation in and out of equilibrium and, by introducing phenomenological 
spin- and charge- fluctuation distribution functions in the computation of the current, obtain the 
differential conductance for large ^. The Hamiltonian to describe the quantum dot system is the 
Anderson impurity Hamiltonian and the current and dot occupation as a function of voltage are 
obtained numerically. We also vary the gate voltage and study the transition from the mixed valence 
to the Kondo regime in the presence of a non-equilibrium current. We conclude with the difficulty 
we encounter in this model and possible way to solve them without resorting to a phenomena 
method. 

PACS numbers: 72.63.Kv, 72.15.Qm, 72.10.Fk 
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I. INTRODUCTION 

The past few years have witnessed a spectacular 
progress in the fabrication and exploration of nano- 
structures giving experimentalists unprecedented control 
over the microscopic parameters governing the physics 
of these systems. Nano-structures, beyond their practi- 
cal applications, display an array of emergent phenom- 
ena stemming from their reduced dimensionality which 
enhances quantum fluctuations and strong correlations. 
Often, experiments are carried out under non-equilibrium 
conditions, with currents passing through the structures. 
The measurements are performed over a wide range of pa- 
rameters, such as temperature and applied bias, allowing 
experimental exploration of the interplay between non- 
equilibrium dynamics and strong correlation physics^"— . 
A canonical example is the non-equilibrium Kondo effect 
observed in a quantum dot attached to two leads held at 
different chemical potentials /i^. The voltage difference 
V = jJ-i — induces a non-equilibrium current I{V) 
through the dot, interfering with and eventually destroy- 
ing the Kondo effect as the voltage is increased. 

The standard theoretical description of the transport 
trough a quantum dot is the two-lead Anderson impurity 
model under a bias voltage. The f- or 2-lead Hamilto- 
nian at zero bias is exactly solvable via Bethe Ansata^!^. 
Using this exact solution as well as NRG calculations for 
example, the thermodynamics of the model have been 
studied in great detail. But the non-equilibrium situa- 
tion, namely when the two leads experience each a dif- 
ferent chemical potential, is much more difficult. This is 
due to the subtle interplay between the non-equilibrium 
aspect of the problem on one hand and the presence of 
strong interactions on the other hand. Technically speak- 
ing it is very non trivial task to find a basis of states that 
diagonalize simultaneously the voltage term and the in- 
teraction term in the Hamiltonian. Nevertheless, a lot of 



efforts have been put forward to study this model but so 
far, only approximate ways of dealing with the voltage 
and/or the interactions have been developed^"^'^ . 

In this paper we develop a phenomenological ap- 
proach to the problem, based on the Scattering Bethe 
Ansatz (SB A), recently developed by P. Mehta and 
N. Andrei (MA)^, a non-perturbative implementation 
of the Keldysh formalism to construct the current- 
carrying, open-system scattering eigenstates for the two- 
lead nonequilibrium Anderson impurity mode. The basic 
idea of the SBA is to construct scattering eigenstates of 
the full Hamiltonian defined directly on the infinite line 
and match the incoming states by two Fermi seas describ- 
ing the initial state of the leads. The non-equilibrium 
steady state transport properties of the system are then 
expressed as expectation values of the current or dot 
occupation operators in these eigenstates. This pro- 
gram has been implemented for the Interacting Reso- 
nance Level Model (IRLM) , a spinless interacting model, 
described in Ref. 7 where the zero temperature results for 
current and dot occupation {hd) for all bias voltages were 
presented. Another exact solution of this model at the 
so-called self-dual point^ by E. Boulat, H. Saleur and 
P. Schmitteckert in Refs. [2^2^ uses conformal field the- 
ory techniques and compares successfully with t-DMRG 
results. 

The main motivation of the present paper is to test 
the very interesting ideas behind the SBA framework on 
a physically more relevant model such as the Anderson 
impurity model and to focus on the phenomenology that 
can be extracted from it. Carrying out the program for 
the non-equilibrium Anderson model we find difficulties 
in the direct application of the SBA approach due to the 
fact that the ground state in the Bethe basis consists of 
bound pairs of quasi-particles, leading to problems in the 
computation of the scattering phase shifts for the quasi- 
particles with complex momenta. This problem is not 
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present in the IRLM when the Bethe momenta are below 
the impurity level and no bound states can be formed. 
We circumvent this difhculty by means of the follow- 
ing argument: The transport property computed in the 
IRLM is related to the single particle phase shift across 
the impurity in the Bethe basis. Based on the same idea 
we develop a phenomenological approach to describe the 
transport property in the Anderson impurity model. We 
identify two types of possible phase shifts across impu- 
rity, which we refer to as " spin- fluctuation" and " charge- 
fluctuation" types to label two phenomenological phase 
shifts akin to the fundamental excitations described in 
the traditional Bethe Ansatz in this model. The phe- 
nomenological Ansatz is checked against exact results on 
the dot occupation in equilibrium and the Friedel sum 
rule21i^, in the linear response regime. Subsequently, 
we discuss our results for the out of equilibrium current, 
conductance and dot occupation. The scaling relations 
for the conductance, predicted from the Fermi liquid pic- 
ture of the problem at strong and weak coupling, are also 
discussed. 

The paper is organized as follows. We start with a 
formal construction of scattering eigenstates in the two- 
lead Anderson impurity model. Then we discuss how 
we impose boundary conditions, which serve as initial 
condition in the time dependent picture, on the elec- 
trons within the leads. Next we shall discuss our re- 
sults for the dot occupation in equilibrium and the con- 
ductance in the linear response regime. Based on the 
checks in equilibrium we then extend our computation 
to the out of equilibrium regime. The difhculty we en- 
counter for complex momenta and the way we handle 
it will also be addressed there. Comparison with an- 
other attempt of exact solution for this model by R. 
M. Konik et al^^^ with the idea of dressed excitations 
above Fermi energy in the Bethe Ansatz picture, first 
considered for the exact conductance of point contact 
device in the FQHE regime^i^, will be discussed. We 
will also comment on the validity and implication of our 
numerical results, among them the charge susceptibil- 
ity, in the out-of-equilibrium regime. Qualitative agree- 
ment between our theory and experimental result is then 
presented. The limit of ?7 — oo is also summarized in 
the last section based on the same phenomenological ap- 
proach. Finally, we summarize our results and conclude 
with some issues on the SBA approach to this model, and 
state how they could be overcome. 



II. THE SCATTERING BETHE ANSATZ 
APPROACH 



A. Scattering state construction 

In this section we apply the SBA approach to con- 
struct the scattering states of the full Hamiltonian. The 



(unfolded) 2-lead Anderson impurity Hamiltonian reads, 

1=1,2'' 

+ W4aiO)d. + dU^AO)) + Ud\d^dld^ (1) 

where summation over the spin indices a is implied. The 
fields ipia- (x) describe chiral, right- moving electrons from 
lead i, U is the on-site Coulomb repulsion between elec- 
trons on the dot, U is the coupling between the dot and 
the lead i, and is the gate voltage. We have set the 
Fermi velocity vp = 

The model's equilibrium properties have been studied 
in great detail via the traditional Thermodynamic Bethe 
Ansatz (TBA)^^. The SBA exploits in a new way the 
integrability of the Anderson Model to construct current- 
carrying scattering eigenstates on the open line. There 
are two main requirements: One is the construction of 
scattering eigenstates with the number of electrons in 
each lead conserved prior to scattering off the impurity. 
Another is the asymptotic boundary condition: that the 
wave function of the incoming electrons, i.e. in the region 
(x <C 0), tend to that of two free Fermi seas far from the 
impurity-. All information about the external bias ap- 
plied to the system is encoded in the boundary condition 
by appropriately choosing the chemical potential of the 
incoming Fermi seas. As in all Bethe- Ansatz construc- 
tions, the full multi-particle wavefunction is constructed 
from single particle eigenstates (now on the infinite open 
line) and the appropriate two-particle S-matrices. We 
first rewrite Eq. ([T]) in the even-odd basis as 



H — He + Ho 



H, = 



dx ipl„{x){-idx)ipeaix) +eddlda 



+ i(C(0)rf. + 4V'e.(0)) + Ud\d^dld^ 



Ho = 



dx ipl„ {x){-idx)ipoa- (x) 



With 



4'eaix) 



^oa(x) = 



Vt'i + t'i 

hlplajx) - tilp2a{x) 



and t = \Jt\ + In what follows we consider the 
case = i2 = for simplicity. The single par- 
ticle solution for even and odd basis is: \e,'pa) = 
fdx{e'P-gpix)i:lix) + epSix)dl)\0) and \o,pa) = 
J dxe'^P^hp{x)^pl^{x)\0), with |0) the vacuum state and 
gp{x), hp{x), Bp independent of spin and given by 

gp{x) = 0{~x) + e'^'-eix) + Sep 0ix)0{-x) , 

hp{x) = e{-x) + 9{x) + Sop e{x)e{-x) , (2) 

_ ^(1 + e'^" + Sep/2) 
2(p-e,) 
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Here 5p = 2 tan ^ ( ) is the usual single particle scat- 
tering phase shift of the electrons off the impurity ob- 

,2 

tained when setting s^p = 0. F = ^ is the width of the 
resonance level. We adopted a symmetric regularization 
scheme 9{±x)S{x) = ^S{x) and imposed \p\ < D, D be- 
ing the bandwidth cut-ofSi. The s{x) = 9{x)9{—x) term 
is a local constant {dxs{x) = 0) in this scheme and it is 
included in the odd channel function to allow the same 
two particle S- matrices, Eq.(l4]), in all channels^^^. The 
9{x)9{—x) term in the even channel wave function is in- 
troduced in order to modify the original (when Sep = 0) 
single particle phase shift across the impurity. The choice 
of Sop and s^p will be addressed later. In the lead ba- 
sis, |z,_pcr), the single-particle scattering eigenstates with 
the incoming particle incident from lead «, can be re- 
stored by taking a proper linear combination of even-odd 
states. For example, \l,pa) = ^(|e,pcr) -|- \o,pa)) — 



jAxi 



(x)|0) is written as 



\l,pa) = j dxe^^^[[9{-x) + \{e^'' + l)9{x)]4,{x) 
+ i(e'^- - mx)i;l{x) + epdUix) + (3) 

with|2,pa) = ^{\e,pcj)~\o,pa)) ^ S Axe'P-al^,{xm 
and s\ {x) related to the 9{x)9{~x) terms by 



!^,/,t 



V2 



V2 



X 9{x)9{-x) 



and 



~ep 



V2 



^op I ■\ 



V2 



X 9{x)9{-x). 



These states have a single incoming particle {x < 0) 
from lead i, that is reflected back into lead i with ampli- 
tude, Rp — {e^^p + l)/2 and transmitted to the opposite 
lead with amplitude Tp = {e'^^" — l)/2. Similar single 
particle states are discussed in Ref. 0. 

The multi-particle Bethe-Ansatz wave-function is con- 
structed by means of the two-particle S-matrix, S(p, fc), 
describing the scattering of two electrons with momenta 
p and k. By choosing Sop — —4 in Eq. ^ (the choice of 
Sep will be discussed in section B and does not affect the 
result here) in the single particle states we can construct 
the same two-particles S-matrix for all combinations in 
even-odd basis (see Appendix. |B| . The two-particles so- 
lution for both particles coming from lead 1 in spin singlet 



state takes the following form 

= J da;idx2^{e*('="^+P"^)Zfep(xi-a;2)al,._^(xi)alp_^(x2)}|0) 

= I J dxidx2A{g{xi, X2)iPI^{xi)^pI^{x2) 

+ h{xi, X2)^l;l^{xi)'il;l^{x2) + j{xi, X2){iIjI^{xi)iPI^{x2) 

-^lix^)4'lAx2))}+ J dxA{e{x){4>lAx)dl~i;l{x)d\) 

+ o{x){i;l^{x)dl - ^Pl{x)d\)) + Amd\dl}\0) 

Here A is the antisymmetrizer. A is an overall normal- 
ization factor and 

g{xi,X2) = Zkp(xi - X2)gk{xi)gp{x2) 
+Zkp{x2 - xi)gkix2)gp{xi) 

j{xi,X2) = Zkp(xi ~ X2)gk{Xl)hp{x2) 

+ Zkp{x2 - Xi)hk{x2)gp{xi) 

h{xi,X2) = Zkp{xi - X2)hk{xi)hp{x2) 
+Zkp{x2 - xi)hk{x2)hp{xi) 
e{x) = Zkp{~x)gp(x)ek + Zkp{x)gk{x)ep 
o[x) = Zkp{-x)hp{x)ek + Zl°{x)hk{x)ep 

m = Zkp{0)ekep 

with Zkp{x) = e~"l"'p9{-x) + e'^'"'9{x) and Zkp{0) = 
ZkpiO). Here tan(</.,p) = ^ 



k+p-U-2ed^'^P^"'- -'^^^^ uc^LL^YKpJ (^k~p){p+k-U-2td)' 

The derivation and more general form of two particles 
case is written in Appendix. [BJ To include spin triplet 

case we denote Zki,kj{xi - xj) = Zk^,kj{xi - Xj)lla'j = 

'Cla',9{xj - Xi) + Sl]l]{ki,kj)9{xi - Xj) where a; is the 
spin index before the scattering and a.^ the spin index 

after the scattering. lala] is the identity matrix. The 
S-matrices must satisfy the Yang-Baxter equations 

sf;°kA:i,fc2)sf;^3(fc^,fc3)s44(fc2,fc3) 

= 8:^(^2, fc3)sfj:53(fc^,fc3)sfi4(fci,fc2) 

for such a construction to be consistent. The two- 
particles S-matrix for this two-lead Anderson model is 
given by 

^_ iB{k}-Bip))■C^+^2UTpii 
^"•"^ " Bik)-Bip)+z2Ur 

with B{k) = k{k - 2erf - C/), P = i(I • I + CT • a) the 
spin exchange operator with and aj representing the 
incoming spin indices. Since the S-matrix is the same 
for all even-odd combinations the S-matrix does not de- 
pend on the lead index i, and the number of electrons in 
a lead, Ni, can change only at the impurity site. This 
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circumstance allows us to construct the fully-interacting 
eigenstates of our Hamiltonian characterized by the in- 
coming quantum numbers, Ni and N2 the numbers of 
incident electrons from lead 1 and 2 respectively. These 
quantum numbers are subsequently determined by the 
chemical potentials fii and /i2. 

To complete the construction of the SBA current- 
carrying, scattering eigenstate, we must still 
choose the " Bethe- Ansatz momenta" {pi}j^^^^ of the 
single particles states to ensure that the incoming parti- 
cles look like two Fermi seas in the region x < 0. This 
requirement translates into a set of " free- field" SBA equa- 
tions for the Bethe- Ansatz momenta-density of the par- 
ticles from the two leads^. The argument is as follows: 
Away from the impurity \i,p(T) reduces to iplai^) with the 
inter-particle S-matrix Eq. Q present. Thus the scatter- 
ing eigenstates describing non-interacting electrons are 
in the Bethe basis rather than in the Fock basis of plane 
waves. The existence of many basis for the free elec- 
tron is due to their linear spectrum which leads to de- 
generacy of the energy eigenvalues. The wave function 
^ipixi+tp2X2^0(^^^ _ 2-2) + se{x2 - xi)]A is an eigenstate 
of the free Hamiltonian for any choice of S with, in par- 
ticular, S = 1 defining the Fock basis and S given in 
Eq. (HI defining the Bethe basis. The Bethe basis is the 
correct "zero order" choice of a basis in the degenerate 
energy space required in order to turn on the interac- 
tions. We proceed to describe the leads (two free Fermi 
seas) in this basis. 

We consider the system at zero temperature and zero 
magnetic field in this paper. To describe the two Fermi 
seas on the leads translates to a set of Bethe Ansatz 
equations whose solution in this case consists of com- 
plex conjugate pairs: p^{X) ~ x{\) ± iy{X) in the A- 



parametrizatioi*i 



33.34.36 



with 



with £d = £d + U/2. Each member of a pair can be ei- 
ther in lead 1 or in lead 2, since the S-matrix is unity in 
the lead space. There are, therefore, two possible con- 
figurations for these bounded pairs. One possible way 
of forming bounded pairs is described by four types of 
complex solutions whose densities we denote aij{X) with 
{ij} — {11, 12,21,22} indicating the incoming electrons 
from lead i and lead j. The other possibility, which is 
perhaps more intuitive in comparing with the free elec- 
tron in the Fock basis, is to include only {ij} = {11, 22}. 
These two types of states give the same results when 
evaluating the expectation value of the dot occupation in 
equilibrium. However when we turn on the bias voltage, 
the results obtained from a 4-bound states description 
show some charge fiuctuations even way below the impu- 
rity level which is not expected from the non-interacting 



{U — > 0) theory (shown in Appendix Thus we shall 
disregard the 4-bound states solution on physical grounds 
and focus on the 2-bound states description in the follow- 
ing discussion. 

To describe in the Bethe basis the two leads as two 
Fermi seas filled up to /ii and fi2, respectively, these den- 
sities must satisfy the SBA equations. 



2(7, (A) = 



1 da;(A) 



dA 



0{X - B,) 

/ dA'if(A-A')^,(A') (6) 



withi^(A) = i^5^^2^. 

The SBA equations are derived from imposing bound- 
ary condition in the free leads (incoming state) region 
and the value of momenta is connected with spin rapid- 
ity A by using the quantum inverse scattering method. 
The Bethe Ansatz equations solved with periodic bound- 
ary conditions at the free lead region with total number 
of particles N {N — Ni + N2 as sum of particle num- 
ber from lead 1 and 2) and the total spin projection S 
(S = Si + S2 ^ N/2 - M with M ^ Mi+ M2 as number 
of down spin particles from lead 1 and 2) are given by 



Xa, + iur 



e J 



Xc - iur 



i = l,2j = l ^ . 



^\ B{k 

^ B{k\)^ X^-iUT 
Xa + iUT 



(7) 



M 

n 



Xa- Xfs + 2iur 
Aa - A/3 - 2iur 



with total energy E = Ei + E2 and Ei — J^j ^\ indicat- 
ing the energy of the electrons within the lead I at zero 
temperature. 

The spectrum of Eq. ^ for one lead case has been an- 
alyzed by N. Kawakami and A. Okiji'^'* where they found 
that the ground state at zero temperature is composed 
of real A; and complex fcj in the thermodynamic limit 
for [/ > 0. The same situation also occurs in the special 
limit where C/ — ^ 00 where P. Schlottmann^^ has done 
also in the one lead case. The proof for two leads ground 
state is similar to the one lead case and is shown explic- 
itly for the finite temperature calculation for the infinite 
U case in Ref . 

As has been mentioned above in the zero temperature 
zero magnetic field ground state all A^ are real (and dis- 
tinct) and fcj form bound state for j = 1,..,2M with 
bound state momenta given by the poles or zeros in the 
S-matrix defined in Eq. Q 

B(k'^{X,)) - A, ± lUV = B{x{X,) ± iy{X^)) -f 7*(A,)(8) 

where 7^ = 0(exp(— L)) and a;(A) and y(A) are shown 
in the Eq. (O. 

Note that the bound state can be formed from four 
possible configurations for B2 < Xa < 00 which we de- 
note bound state from lead i and lead j quasi momenta 
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denoted as . The bound state between i?i < Aq < B2 
can only be formed by quasi momenta both coming from 
lead 1. As already mentioned the four bound state dis- 
tribution does not give physically sensible results for the 
charge susceptibility as shown in Appendix A and there- 
fore we will limit our discussion to two types of bound 
state distribution here. Below we surpass the index of 
lead in A and put back the index dependence in the end 
for simplification. Inserting Eq. ([8]) into Eq. (O we get 



gjfc+L ^ TT K- Xp + 2iUT 
jJi A„-A^+7+ 

ifc-L ^ TT -^a - + la 

J-i K-Xb- 2iUT 

13=1 ' 
-pj Aff - Aa + 7^ _ ^ 

f}^ Xfi-K + ip 



(9) 
(10) 
(11) 



Thus for L — > 00 from multiplication of Eq. Q and 
Eq. (dni) we have 



i _ -pr Aa - A/3 + 2iUT 
2iUT 



n 



(12) 



+ 7r (13) 



Taking the logarithm of Eq. ([T2|) we have: 
2nJa = -2x{K)L J2 (2^2 (^^1^ 



with Onix) = tan~^(2x/n) and {Ja} a set of integer num- 
bers. We can extend the definition of Ja to include inte- 
gers or half integers and rewrite Eq. (|13p as 



TT 



-x[\o 



Xg — Xfz 

2UT 



(14) 



Now let us put back the dependence in lead indices. 
Starting from Eq. ([14]) it can be shown that there is 
one-to-one correspondence between the Aq.'s and the Jq.'s 
and that all Ac's have to be different. Thus the set 
of rapidities {A^}, characterizing an eigenstate of the 
Hamiltonian, is uniquely determined by one specific set 
of {Jq}. For instance, the ground state of the Hamilto- 
nian Hq in the presence of a bias voltage is simply ob- 
tained by packing two "Fermi seas" of non-consecutive 
integers (Pauli principle in lead space) up to certain 
"Fermi points" (see Fig. [T] ) corresponding to the Bi 
and B2 in the continuum limit. For notational simpli- 
fication we relabel {ij} = {11,22} as {/} = {1,2}. Now 
defining Y.^,<>^i) = Tit ^ Ei'^^'HK) and using 
dxOn{x) = i_|_(2a:/n)^ "^^^ Write Eq. (jl4p in the contin- 
uum limit (by taking L — >■ 00 and differentiate Eq. 
with respect to A). Doing so we shall distinguish two 
different domains: 

For i?2 < A < 00 the particles are fully packed and 
states are labeled by a different lead index /. In this 



lead 1 



lead 2 



00 



Bi 



00 



B2 



FIG. 1: Sketch of the configuration of Bethe momenta corre- 
sponding to the ground state of Ho with an additional bias 
voltage i.e. two Fermi seas at different chemical potential. 



domain, the SBA equations in the continuum limit takes 
the form 



E-^'nA) = -i^- / dA'i^(A-AV^)(A') 



1=1 



B2 



dX' K{X- y)a^^HX') . (15) 



Bi 



For Bi < A < i?2 we can see from Fig. [T] that the lead 2 
states are unoccupied. We shall introduce a distribution 
of holes for the lead 2 that we win denote 5-(2)(A). The 
continuum SBA equations in this regime are given by 



aW(A)+a(2)(A) = -i^- 
TT dA 



dX' K{X-X')a^^\X') 



dX' K{X-X')a(^\X') (16) 



Since a^'^\X) obeys the same equation as (t(2)(A) as may 
be seen from subtracting Eq. (|T6| and Eq. (|15p we can 
combine Eq. ([T5|) and Eq. p^ together to get 

2a(A) = -i^^-2 r dX' K{X-X')a{X') 
dA Js, 



dA'if(A- A')ct(A') (17) 



with Bi < A < 00 for lead 1 and ^2 < A < 00 for lead 
2 Bethe momenta density distributions. Each density is 
defined on a domain extending from Bi to the cutoff D 
- to be sent to infinity. The B^ play the role of chemical 
potentials for the Bethe- Ansatz momenta and are deter- 
mined from the physical chemical potentials of the two 
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leads, fj,i, by minimizing the charge free energy, 



dX (x(A) - iJLi)(7i{\) 



(18) 

with (Ti the lead 1 particle density and 1J2 the lead 2 
particle density. Note that ai and 12 obeys the same 
integral equation Eq. ^ with different boundary (cri(A) 
with A C (i3i,oo) and cr2(A) with A C (i?2,oo)). Solving 
the SBA equations subject to the minimization of the 
charge free energy fully determines the current-carrying 
eigenstate, and allows for calculation of physical 

quantities by evaluating expectation value of the corre- 
sponding operators. In the following we shall discuss our 
results from equilibrium cases to non-equilibrium ones, 
starting with the expression for various expectation value 
of physical quantities. 



object is because Bethe wave- functions are orthogonal to 
each other for different pairs of Bethe momenta under the 
condition that the size of the leads L taken to infinity. 

Here we address the different choice of Sep (with Sop = 
—4 fixed to have the same S-matrix in all channels) which 
gives rise to different forms of v{X). We shall first discuss 
the natural" choice Sep = (i.e. absence of 0{—x)9{x) 
terms) and show it reproduces the exact result for the 
dot occupation in equilibrium. While in checking the 
steady state condition, i.e. d{nd)/dt — 0, for out-of- 



equilibrium situation, the choice of Sg 



fails. To 



remedy this issue we propose Sep (i.e introducing 
counter-intuitive 0{—x)9{x) terms) schemes to circum- 
vent this difhculty. We check this proposed phenomena- 
logical scheme in equilibrium against the exact dot occu- 
pation obtained in Sep = case in the second part of the 
discussion as a benchmark for our approach. First let us 



discuss the result for Se 



0: 



(1) Se 



0: We choose Se 



B. Expectation value of current and dot occupation j^^^ Anderson impurity model. Denote v{X) 



as in the case of the 1- 
(A) 



,SBA 



For /ii = /i2 all Bi are equal to some equilibrium 
boundary B fixed by the choice of /i^ . The dot occupation 
is given by the expectation value '^^^{'^ , iii\d\,da\^ , iJ-i) ■ 
Taking the limit L — 00 (L being the size of the 
lead) one can express n^t as an integral over the den- 
sity of A and the corresponding matrix element ^{X) ~ 

^'"''/'.'A?4-^fltAr-mr'''^ taken to order f Here 

(p+(A)p (A)|p+(A)p (A)) L 

the state \p'^{\)p~ (A)) denotes a pair (or bound states) of 
quasi-particles with complex momenta given by Eq. ([S]). 
The reason why Ud is governed solely by one-bound state 
matrix element instead of a complicated many-particle 



in this choice. The dot occupation expectation value in 
equilibrium is given by 



nd = 



^i2\J2adU^\'^,^^l= 1^2) 



M2|*,Ml 

= 2 



dXa{X)iy^^^{X) (19) 



where the factor 2 in front of the integral accounts for 
the spin degeneracy. The matrix element of the operator 
dl^da in the SBA state is given by 



J 



SBA 



(A) 



2r 



162;(A)r2 



xHX)+yl{X) [xHX)+yUX)WW+ylW] \2x{X)-U 



r 



where we introduced, for simplified notations, the func- 
tions x(A) = x{X) — ed and y±{X) — y{X) ± F. 

Eq. (jl9p can be proved to be exact by comparing it 
with the traditional Bethe Ansatz (TBA) result. In the 
latter, Ud is computed as the integral of the impurity 
density. This observation that the SBA and TBA results 
for Ud agree in equilibrium shows the connection between 
the dot occupation and the dressed phase shift across 
the impurity. The dressed phase shift mentioned here is 
equivalent to the impurity density as can be seen in the 
Eq. (|C1I) in Appendix [Cj The proof of the equivalence 
between TBA and SBA in equilibrium is also given in 
Appendix [C] 

To describe the out-of-equilibrium state we first check 
if the steady state condition = (or equivalently. 



d{Ni_+N2) 
dt 



0) is satisfied in this basis. As mentioned ear- 



lier these scattering states are formed by bounded quasi- 
particles with complex momenta and therefore the single 
particle phase across the impurity is not well defined in 
the sense that |e'*p* I 7^ 1- This problem begins to sur- 
face as we set out to evaluate transport expectation value 
and renders 



dt 



B22 



Bii 



with 



A(A) 



dAcr6(A)A(A) ^ 



(20) 



x^X)+e{XmHX)+yliX)] 



Thus it appears that using this basis the steady state 
condition is not observed. This problem does not appear 
when the momenta are real as in the IRLM case^. 
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(2) Sep 0: To remedy this problem we redefine the 
single particle phase shifts across the impurity, in anal- 
ogy to the results for the IRLM^, through the choice 
of nonzero Sep in Eq.Q. With a suitable choice of 
Sep we may restore a well defined single particle phase 

[e^'^p*! — 1 with Sp± denoting this new phase. The way 
we judge whether we make the correct choice for the new 
phases 5p± is to compare the dot occupation Ud in equi- 
librium before and after the redefined phase. The ex- 
plicit form of Sep and phase Sp± will be motivated below 
but first we shall show that a single redefined phase is 
not sufficient to satisfy the constraint of dot occupation 
comparison. 

Again the choice of new phases is constrained by the 
requirement that we shall obtain the same result for 
{J2cr^tda} as given by z/'^^^(A) in equilibrium. Based 
on this constraint it can be shown explicitly that a single 

well defined phase (in the sense of \e'^p^ I = 1) is not suffi- 
cient to reproduce the equilibrium v^^^lX) as following: 
The new dot amplitude ep+ and Cp- have to satisfy 



lep+l' + IS-P 



4r 



e„+ e. 



|2 



4r2 



As both |ep+p and |ep- p are positive we see that a sin- 
gle redefined phase cannot satisfy the above constraints 
simultaneously. Therefore we have to choose at least two 
sets of redefined phases Sp± (with i = s,h denoting spin- 
fluctuation or charge-fluctuation to be addressed later) 
and, along with them, some distribution functions /* to 
set the weight for these phases. 

To motivate the idea of searching the correct phase 
shifts we shall come back to the derivation of dot occupa- 
tion in traditional Bethe Ansatz (TBA) picture. In TBA 
the total energy of the system is described by energy of 
the leads electrons and energy shifts from the impurity. 



(21) 



Based on Feynman-Hellman theorem, which is applicable 
in equilibrium (closed) system, we have 



ded L ded L 

3 3 



ded 



The result for Eq. (|22l) agrees with those obtained from 
Eq. (jC2l) and can be viewed as a third approach to ob- 
tain the expectation value of the dot occupation. The 
key observation here is that this quantity is related to 
the bare phase shift Sp+ +Sp- and therefore the redefined 
phases must be proportional to this quantity. Among 
them there are two likely candidates with redefined phase 
shift given by Sp+ + Sp- , describing the tunneling of a 

bounded pair, and -^^i^— , describing the tunneling of a 



single quasi-particle. In a sense this is the echo for the ele- 
mentary excitations above the Fermi surface in the Bethe 
basis characterized by N. Kawakami and A. Okiji'"^ as 
charge-fluctuation excitation, which describes bounded 
pair quasi-particles excitation, and spin-fluctuation exci- 
tation, which describes one quasi-particle excitation. An- 
other similar picture is the spin-fluctuation and charge- 
fluctuation two fluids picture proposed by D. Lee et al^ 
albeit in a different context. We identify the phase de- 
fined by 

? ^ Sp+ + Sp- _ ~ 
dp- = dp+ = = dp 



(with Sg 



„,=l(z(p±-e,)-r)(e^(--^)-l)) 



as spin-fiuctuation phase shift and 



Sp+ 



(with s 



ep^ 



^ep± 



^(»(p±-erf)-r)(e'(V+V)-l)) 
as charge-fluctuation phase shift. 

The out-of-equilibrium current is evaluated by the ex- 
pectation value of current operator / with (/) defined 

by 



{!) = 



mA0^)~4A0^))d.~h.c.)) (23) 



in the state Notice that tpj^^O''^) = 

lime-^Q{tpl^{—e) +i/'L(+^))/2 is introduced in transport 
related quantity to be consistent with our regularization 
scheme which introduces another local discontinuity in 
odd channel at impurity site. 

From Eq. and the expression for the phases 5p and 

Sp we have the expression for current as 

^(M1,M2) = (*,Ail,Ai2|/|*,Ml,M2) 



2e 



dA abiX)ifs{X)r{X) + fh{X)J'\X)) (24) 



Bi 



The corresponding spin-fiuctuation and charge- 
fiuctuation matrix element of the current operator based 
on the spirit of Landauer transport, denoted as J'^(A) and 

J''(A) with J"(A) = |Tp(A)|2 = (a = {s, /i}) de- 

pending on redefined phase shift Sp only, are given by 

. , , sgn(x(A))(.^(A) + „^(A)-n 

V(i2(A) + y2(A) - r2)2 + 4r2£2(A) 
2r2f2(A) 

(x2(A) + r2)2 - 2y2(A)(r2 - 52(A)) + y4(A) ' 

(26) 



J\X) = 



Here sgn(a;) = is the sign function. It is introduced 
in order to pick up the correct branch when taking the 
square root in denominator of Eq. (1251) . This way we en- 
sure that J* (A) has the proper limit when U is sent to 
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infinity (cf Section III) . Other than the motivations men- 
tioned above for identifying spin and charge fluctuation 
phase shifts the functional forms of J" (A) and J'^iX) as a 
function of bare energy a;(A) can also be used to identify 
these two type of phase shifts (See Fig. [12] in Section III 
for infinite U Anderson model, the finite U is similar). 

Next we shall choose the appropriate weight for each 
type of phase shift. So far we have not yet been able 
to deduce the form of these weight functions /s(A) and 
//j(A) and we introduce them phenomenologically. Let 
us define phenomenological spin-fluctuation and charge- 
fiuctuation weight functions as 



and 



hieiX)) 



Dh{e{X)) 
Ds{e{X))+Dh{e{X)) 



(27) 



(28) 



Here Ds{e{X)) is the spin-fiuctuation density of state, 
Dh{e{X)) is the charge-fluctuation density of state as de- 
fined in Ref. 113, and £(A) is the corresponding dressed 
energy i.e. the energy required to produce these spin- 
and charge-fluctuation excitations above the Fermi level. 
Here dressed energy refers to the sum of the bare en- 
ergy of adding/removing one bound state, as in charge 
fluctuation, or single quasi particle, as in spin fluctuation, 
and the energy shift from other quasi particles due to this 
change. The equation that solves a single quasi-particle's 
dressed energy £(A) reads^ 



e(A) = (xiX) - m) 



dX' K{X- X')e{X') , 



(29) 



We wish to compare at this point our approach to the 
one taken by Konik et a P^'^° . The authors' Landauer 
approach is based on an ensemble of renormalized exci- 
tations, the holons and spinons, and the conductance is 
expressed in terms of their phase shift crossing the impu- 
rity. However, the leads are built of bare electrons and 
thus one faces the difficult problem of how to construct 
a bare electron out of renormalized excitations in order 
to be able to impose the voltage boundary condition. 
The basic approximation adopted, electron antiholon 
+ spinon, is valid only when the electron is close to the 
Fermi surface (see N. Andrei^), and therefore the ap- 
proach is trustworthy only for very small voltages. Nev- 
ertheless, the dressed excitations framework seems to give 



at least qualitatively good results when another energy 
scale (such as the temperature or an external field) is 
turned on^. In contrast we construct the eigenstates of 
the Hamiltonian directly in terms of the bare electron 
field and can therefore impose the asymptotic boundary 
condition that the wave function tend to a product of 
two free Fermi seas composed of bare electrons. While 
we do not have a mathematically rigorous derivation of 
the weight functions we introduced, the validity of the 
scattering formalism is not restricted to any energy win- 
dow other than energy cutoff. 

C. Results for equilibrium and linear response 

In the numerical computation, for the practical pur- 
pose, we assumed Kondo limit (U = —2ed, ^ ^ 1) form 
of the spin-fluctuation and charge-fluctuation distribu- 
tions, i.e. 



^e\X)+Tl 



(30) 



and 



Dh{e{X)) 



1 



V2t7r (e(A) + e<i)2 + r2 
with Tj^ being the Kondo scale derived in Ref. |40i as 



(31) 



Tk = e^ 2t.r 



(32) 



As we use the Kondo limit in our expression for the 
spin-fluctuation and charge-fluctuation distributions, we 
expect our phenomenological approach works better for 
large U/T. We also take e(A) ~ x{B) — x{X) for numer- 
ical convenience with B denoting the Bethe momenta 
boundary given by /ii = /i2 =0. The dot occupation 
(Tlia '^ct^ct) evaluated by these new phases is given by 



(^4d,) =2( j^^ dXab{X){u'{X)fs{X)+v\X)h{X)) 



+ dAab(A)(//^(A)/,(A) + i.'*(A)A(A)) I (33) 



with z^'*(A) and i^''(A) given as 
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1 - 



v/(i2(A) + y^{X) - r2)2 + 4r2i2(A) 



^ 'T \ ^(.i2(A) + y^{\) - r2)2 + 4r2i2(A) 



i(A) 



2i(A) - [/ 



(34) 



2ri2(A) 



_ (.i2(A) + r2)2 _ 2y2(A)(r2 - i2(A)) + y4(A) 

1 



36?;(A)ri2(A) 



(£2(A) + r2)2 - 2y2(A)(r2 - i2(A)) + /(a) \^2i(A) - U 



(35) 



respectively. We may check whether this choice of phe- 
nomenological distribution functions satisfy the condition 
in equilibrium that 



(V4d,)=4/ dAafc(A)^.^^^(A) 

JB 



= 4(^y^ dAab(A)(z.^(A)/,(A) + i.''(A)A(A))j . (36) 

We can see from the Top of Fig. [2] that the compari- 
son between the phenomenological and the exact result 
for the dot occupation in equilibrium is good deep into 
the Kondo regime (e^ ~ ^""^^ away from it 

(ed ^ 0) but is worse when we are in mixed valence 
region (e^ ~ 0). This discrepancy, due in part to the ap- 
proximations we made for Ds(e) and Dh{e), may go away 
if we took more realistic form of Ds{e{X)) and Dh{e{X)) 
also in mixed valence regime as suggested in Fig. [2] How- 
ever the numerical procedure is much more complicated 
there. We confine ourself to this simpler limit in our 
phenomenological approach. 

Another check on our result in equilibrium is to find the 
linear response conductance through our formulation and 
compare with the exact linear result given by the Friedel 
sum rule2'^'28. The Friedel sum rule, which relates the 
equilibrium dot occupation to the phase shift experienced 
by electrons crossing the dot, is related to zero voltage 
conductance by ^|v=o = 2sin^(7r(nd) /2). The zero bias 
conductance in our construction can be analyzed easily'^^ 
by noting that at low- voltage eV ^ ^2 — ~ 

N2) = 47r at{X)dX. By taking B2 ~ Bi = B in the 
expression for the current across the impurity Eq. (j24p 
we get the zero bias conductance expressed as 



dl_ 

dV 



v=o 



[MB)J%B) + MB)J\B)] (37) 



Here B = B{p, e^, F, U) is determined by Hi = H2 = 0. 
The comparison between Friedel sum rule (FSR) re- 
sult and the conductance given by Eq. ([57)) (denoted as 



(pSBA)) is shown at the Bottom of Fig. [5] It displays 
the consequence of the equilibrium Kondo effect in the 
quantum dot set up: due to the formation of the Kondo 
peak attached to the Fermi level the Coulomb blockade 
is lifted and a unitary conductance is reached for a range 
of gate voltages around —U/2. Again we see that the 
comparison is good for large U/T but poorer in mixed 
valence regime for smaller U/T, which is consistent with 
the observation we made when evaluating {fid) as shown 
in top figure of Fig. [51 Having checked our results in 
equilibrium we shall go on to compute the current and 
the dot occupation in the out-of-equilibrium regime. 



D. Results Out-Of-Equilibrium 

Now let us begin to investigate the current and dot oc- 
cupation change as we turn on the voltage. We start with 
the discussion on current vs voltage for various regime. 
The current vs voltage is plotted in the inset of figure of 
Fig.Ofor different values of U and at the symmetric point 
ed = —U /2. Note that we use an asymmetric bias voltage 
when solving numerically the integral equations originat- 
ing from Eq. ([5]) with constraint of minimizing the charge 
free energy Eq. (fT5)) : Namely we fix /ii ~ (around 
lO^'^ — 10~^) and lower ^2- Therefore, a direct confronta- 
tion between the results obtained from real-time simula- 
tions of the Anderson model out-of-equilibriu m^^'2°i23 jg 
difficult but the main features of our calculation match 
the predicted results: a linear behavior of the I-V charac- 
teristics at low- voltage, the slope being obtained from the 
FSR (2 in units of /h at the symmetric point), and a 
non-monotonic behavior at higher voltage, the so-called 
non-linear regime. In particular, our calculations show 
clearly that the current will decrease as U/T is increased 
which is in agreement with other numerical approaches 
(e.g. cf Fig. 2 of Ref. [23 for a comparison). 

The plots of the differential conductance vs source 
drain voltage for different dot levels, e^, tunneling 
strengths F and interaction strengths U are shown in 
Fig. Oand Fig.O Two major features emerge from these 
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defined as 











\ VT 


r=0.5(TBA) - 




o r=0.5(pSBA) . 




1 — U.ZO^ 1 Of\) 




V 1 — U.^O^pOOAJ 




r=0.1(TBA) 




t> r=0.1(pSBA) 











2ur 



(38) 



with ci = 0.002. Tlie energy scale was ex- 
tracted from tlie numerics by requiring tliat the function 
dI/dV{V/T^) decreases to half its maximal value when 
V o^T*. The expression for T* as given by Eq. §^ dif- 
fers from the thermodynamic Tfe as defined in Eq. (j32p . 
The difference of prefactor in the exponential is certainly 
related to the unusual choice of regularization scheme in 
the SBA^. The other possible implication for this dif- 
ferent formulation for the Kondo scale is also addressed 
later when we discuss the experiment done by L. Kouwen- 
hoven et al^. 
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FIG. 2: Top: (rid) as a function of ed from the exact result 
(dotted line) and from Eq. ((36]) (solid line). Bottom: The 
differential conductance in the linear-response regime, as a 
function of ed from the phenomenological Scattering Bethe 
Ansatz (pSBA) and exact linear response conductance from 
Friedel sum rule (FSR) for F = 0.5, 0.25, 0.1, and U = 8. 




V/r 



FIG. 3: dl/dV vs V/T for F = 1, td = -U/2, and various 
U. Inset: Steady state current vs voltage curves for F = 1, 
ed = —U /2, and various U. Dashed line is a line with constant 
conductance =j- plotted for comparison. 



plots: 1) A narrow peak around zero bias reaching max- 
imal value of 2e^/h (the unitary limit) for values of the 
gate voltage close to the symmetric point (e^ ~ —U/2). 
2) A broader peak developing at finite bias. The first 
peak is a non-perturbative effect identified as the many 
body Kondo peak, characteristic of strong spin fluctua- 
tions in the system. But the broad peak is due to renor- 
malized charge fluctuations around the impurity level. 
Notice the two features merge as the gate voltage, ed 
is raised from the Kondo regime, ed = —U/2, to the 
mixed valence regime, ed = 0, with the Kondo effect dis- 
appearing. As a function of the bias the various curves 
describing the Kondo peak for different values of the pa- 
rameters can be collapsed onto a single universal function 
dl/dV = dl/dViV/T*) as shown in Fig. H Here T* is 



The small voltage behavior for differential conductance 
in symmetric case, i.e. ed — — -j, is expected to beiiii^ 

and allows us to identify the constant ay from the 
quadratic deviation from 2e^//i. The quadratic fit of the 
universal curve around y ~ 0, as shown in Fig. 21 gives 
q;v ^ 1. It is also expected ior <^ V <t: ^ that the 
tail of the peak decays logarithmically — as 

dl 2e2 1 

k 

The latter behavior is observed (see inset of Fig.|4]) in the 
regime ^ 1 for 10^ < ^ < 10^ with the logarithmic 
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0.5 



0.0 




V/T. 



FIG. 4: Top: Zoomed in picture of the differential conduc- 
tance vs voltage nearby zero voltage. Inset shows the univer- 
sality in conductance vs voltage scaled by when < 1. 

k 

The quadratic behavior occurs for ^ < 0.5 as indicated by 
the fitted curve. Bottom: Differential conductance vs voltage 
scaled by nearby the Kondo peak structure. Inset shows 



the logarithmic behavior when 
data sets. 



> 1. r = 0.5 for aU these 



function given by 

dl _ 
dV ~ 



C2 



with the parameter C2 = 0.055. Here /(^) is simply 
a constant (in V) shift. As suggested from the bottom 
plot of Fig. 0] (see also Fig. [T3] for the infinite U case) the 
charge fluctuation side peak does not fall into the same 
scaling relation but the strong correlations shift the cen- 
ter of the side peak closer toV = (see Fig.[3]and Fig.[5|). 
In other words the position of the resonance in the dl / dV 
curve naively expected around V — \ed\ is renormalizedift 
by the presence of interactions. In the inset of Fig. [5] we 
show the logarithm of the voltage obtained at half width 
half maximum (HWHM) of the zero voltage peak and 



FIG. 5: dl/dV vs V/4r for i7 = 8, F = 0.25 and various 
ed from Kondo (td = —4) to mixed valence regime (td ^ 0). 
Inset: Comparison of In(r^) — In(ci) and ln{VHWHM) as a 
function of impurity level ed- Here Vhwhm is the voltage 
difference estimated at half value of differential conductance 
at zero voltage. The constant shift — In(ci) is chosen to give 
the best fit in the data away from 



u_ 

2 ■ 



compare it with 



lnT,t^--(--i^i+^%lnfcr-^^ 



2UT 



(after subtracting the constant Inci). What is impor- 
tant and universal is that both quantities (In Vhwhm 
and In ) exhibit a quadratic behavior in the gate volt- 
age e^. Similar results had been found experimentally 
by L. Kouwenhoven et al^ when they compare the full 
width half maximum of dl/dV (from which they obtain 
a Kondo scale Tk^ at finite voltage) with the tempera- 
ture dependence of the linear response differential con- 
ductance (from which another Kondo scale T^^ is ex- 
tracted). It is suggested from our numerical results that 
both InTfcj (in analogy with our Tk) and InT/j^ (which 
is our T^) follows similar quadratic behavior in ed but 
differ in their curvatures by a factor of tt. In Ref. |5| the 
curvatures of the quadratic behavior differ by a factor of 
around 2 (see Fig.SB in Ref. which is attributed to 
dephasing of spin fluctuations at finite voltage. 

Notice that in all the numerical data shown for cur- 
rent vs voltage we have chosen > 8 to explore the 
scaling relation in the Kondo regime. Another reason is 
that our phenomenological distribution functions intro- 
duced to control the relative weight for spin- and charge- 
fluctuation contributions work is much better in the large 
Y regime (cf. Fig. [5]). 

Next let us study the change in the dot occupation as 
a function of the voltage. The extension of the com- 
putation of the dot occupation out of equilibrium is 
straightforward. Suppose we find the correct distribu- 
tion functions fs{X) and /;i(A) then we have iy^^^{X) = 
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FIG. 6: (n'd) vs V/T for different U with ed = - f and T = 1 
case. Inset: The corresponding nonequilibrium charge sus- 
ceptibility. A small peak shows up nearby V = Q for all these 
curves. 



FIG. 8: for various fixed voltages as a function of ea 

for r = 0.25, U — d>. Inset shows {hd) vs td for various fixed 
voltage. 




FIG. 7: vs V/4r for F = 0.25, [7 = 8, and various 

from Kondo to mixed valence regime. We see that the small 
peak nearby V — Q only appears when — >■ Inset: The 
corresponding (fid) vs V/AY. 



v'{\)fs{\)+v^{\)fh{X). Under this assumption v^^^{\) 
retains its form in and out-of-equilibrium and the general 
expression for {hd) is 



(39) 



dA c7f,(A)i/^^^(A) + / dA (7b(A)jy*«^(A) 



As the form for ^^'^^{X) is proved to be exact in equilib- 
rium, we shall regard Eq. ([5^ as an ea;aci result for {hd) 
even out of equilibrium and valid in all different range of 
U, ed, r under the assumption that the integrand does 
not change its form for in and out of equilibrium, which 
is the case for general results of SBA. In the numeri- 



cal results shown hereafter we shall use this expression. 
Eg. pop , for matrix element of dot occupation rather than 
Eq. p6|) . We adopt the same voltage drive scheme by fix- 
ing fix and lowering fi2. 

By using this result we do not need to confine ourself 
for large jr- The case for different ^ with €d — ~y 
and for U = 8,r = 0.25 with different are shown in 
Fig. |6] and Fig. [T] The main features of these plots are a 
relatively slow decrease of the dot occupation at low volt- 
age followed by an abrupt drop of {ud). The decrease of 
{ud) takes place within a range of voltage of the order 
of r. Then as we increase the voltage further another 
plateau develops. Note that, as expected, the bigger U is 
the higher the voltage needed to drive the system out of 
the {ud) — 1 plateau. In a sense the charge fluctuations 
are strongly frozen at large U and it costs more energy 
to excite them. The voltage where the abrupt drop in 
{ud) occurs corresponds to the energy scale at which the 
"charge fluctuation peak" was observed in the conduc- 
tance plots. This can be seen by comparing the position 
of the broader peak in Fig. [5] with that of the abrupt dot 
occupation drop in Fig. [71 

Similar to the differential conductance we may define 
the nonequilibrium charge susceptibility as 



Xciv)U 



d{hd) 
dV 



that we obtain by taking a numerical derivative of the dot 
occupation data with respect to the voltage. In the case 
of J7 = — ed/2 there are two features as can be seen from 
the inset of Fig. [S] and main figure of Fig. [71 Nearby V ~ 
we see a first small peak arising with width and height 
decreasing with increasing ^. We identify this peak as 
a small remnant of the charge fluctuations in the Kondo 
regime. This statement is confirmed by noticing that this 
peak goes away as jr increases, vanishing when ?7 — > oo 
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as shown in Section III where the infinite U Anderson 
model is discussed. The second peak is located at the 
same voltage as the charge fluctuation peak observed in 
the conductance plots and is therefore associated to the 
response of the renormalized impurity level to the charge 
susceptibility. This can be seen when comparing Fig. [5] 
and Fig. [71 

Another interesting quantity, the usual charge suscep- 
tibility, defined by Xc{(-d)\v — — ^glf ■, can also be quali- 
tatively described. In Fig.|S]we plot —^^j- as a function 
of td as we only have a few points in fixed td for finite 
voltage. Notice that Xc{^d)\v tends to be an universal 
curve in large voltage, indicating charge on the dot re- 
mains at some constant value in the steady state with 
large voltage. This constant value at large voltage, as 
pointed out by C. J. Bolech, is around 0.65 for ed — — ^; 
case. In preparing this article we noticed that a simi- 
lar computation, adopting the same asymmetric voltage 
drive protocol as we have here, is carried out by R. V. 
Roermund et al2i for the dot occupation out of equilib- 
rium by using equation of motion method. We do get a 
similar value for the dot occupation at large voltage. This 
value is different from the dot occupation value Ud — 0.5 
at large voltage when the interaction U is turned off as 
shown in Fig. [T31 This difference might have to do with 
the 0.7 structure observed in quantum point contact* in 
high temperature (temperature is high compared with 
the Kondo scale but still small compared with phonon 
modes or electronic level) and zero magnetic field as the 
linear response conductance given by Ud = 0.65 by using 
Friedel sum rule is around 0.73. In a sense the voltage 
seems to play a similar role to the temperature on the 
way it influences the dot occupation. Further connection 
between these two behaviors could be clarifled by com- 
puting the decoherence factor as in Ref. [2l|. This deco- 
herence factor is related to the dot correlation function 
out of equilibrium which can be computed in three-lead 
setup^ by using our approach. 

E. Comparison with other theoretical and 
experimental results 

In most of the other theoretical approachesi^ii^"— i^^i^ 
the symmetric voltage drive (/ii = — /J.2) is usually as- 
sumed to preserve particle-hole symmetry in symmetric 
case {ed = —%)■ It is thus difficult for us to make any 
definite comparison with other theoretical results. The 
qualitative feature, as shown by the black curves in Fig. [9] 
done by D. Matsumoto^S by using perturbation expan- 
sion in U at strong coupling fixed point, is similar to our 
results in the sense that the height of the charge fluctua- 
tion side peak and width are almost the same. The major 
differences are in the shape of Kondo peak and the posi- 
tion of the charge fluctuation side peak. A clear signature 
of renormalized dot level Ed as hinted in renormalization 
computational^ is clearly seen in our result. The shape 
of Kondo resonance nearby zero voltage deviates from its 



quadratic behavior expected from Fermi liquid picture at 
smaller voltage in our case as is expected for asymmetric 
voltage driv o^^i^^ . 

1.0 j — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — I 
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FIG. 9: Comparison of our theory with perturbation expan- 
sion in U done by D. Matsumoto on dl /dV (y-axis in unit of 
2e^/h) vs V/U (x-axis). Our data (Blue, purple, and brown 
lines correspond to ^ = 0.13,0.083,0.063 respectively. A 
shown in inset is F in our notation. EQ in the inset is con- 
ductance computed by equilibrium density of state which is 
not relevant to our discussion here.) is shown as the main 
figure and Fig. 8 in Ref [2^ is shown in the inset. In Ref [2^ 
the voltage is driven symmetrically, i.e. fii = —^2, rendering 
the factor of two difi'erence in the voltage (i.e. = 0.5 in 
our case corresponds to ^ = 1 in the inset, e = 1 in our 
convention.) in comparing our result with that in Ref |2^. 
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FIG. 10: Comparison of theory with experiment of dl /dV 
(y-axis in unit of /h) vs V (x-axis in unit of mV). Inset 
is the original data graph published in Ref. ^ The red dots 
are given by our theory for = 8 with voltage rescaled to 
fit with original data in unit of mV . The value of differential 
conductance (experiment data in black line) is rescaled from 
(0.6,1.3) to (0,2) in unit of ^. 

We can also compare our results with experiments. As 
shown in the inset of Fig. [10] is the ^ vs measured 
in Co ion transistor by J. Park et al.^. We rescaled the 
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differential conductance and superimposed our numeri- 
cal results on the data graph. The measurement was 
done by using an asymmetric drive of the voltage (by 
keeping /ii = and changing /i2 to be larger or smaller 
than zero) and thus there is an asymmetry in the differ- 
ential conductance as a function of voltage as illustrated 
in the data curve. In our numerics we only compute the 
scenario for /ii = and lowering /i2 (only for V > re- 
gion of Fig. [To]). The V < region is plotted by just 
a reflection with respect to the ^ = axis which illus- 
trates the case of /i2 = and lowering ^i. To compare 
with the correct voltage setup on the V < side as in 
experiment will involve computations within a different 
parametrization for bare the Bethe momenta which is be- 
yond our current scope. The comparison on the ^ > 
region shows good agreement between our theory and ex- 
perimental result. The discrepancy on the width of the 
charge fluctuation side peak could be due to the vibron 
mode"**. To describe these type of transistors we shall 
start with the Anderson-Holstein Hamiltonian. We are 
currently exploring the possibility of solving this model 
by the Bethe Ansatz approach. 



III. INFINITE U ANDERSON MODEL 



In the limit of jr 



oo the finite U two-lead Ander- 
son impurity Hamiltonian becomes the two-lead infinite 
U Anderson model. The latter model is closely related, 
via the Schrieffer- Wolff transformation^, to the notori- 
ous Kondo model, a model of spin coupled to a Fermi 
liquid bath. The reason for that is simple: since U ^ oo 
the charge fiuctuations are essentially frozen out and only 
the spin fiuctuations dominate the low-energy physics. 
The Hamiltonian is given by 

i=l,2 

+ U{^l{0)bU, + dlbij,^{Q)) (40) 

Here the bosonic operator b is introduced to conserve 
b^b + J2cr d-l-da = 1 and by applying the slave boson tech- 
nique we project out the phase space of double occupancy 
occurring in finite U case. The corresponding Bethe mo- 
menta distribution function for the infinite U Anderson 
model is given by 



2ct(A) = -- I dA'K{A - A')ct(A') 

TT 



Bi 



dA'K{A - A')ct(A') (41) 



with K{A) ^ 1 (2r)^+fA-A0^ ■ 

Eq. (|4ip can be derived directly following the proce- 
dures in the finite U Anderson model. It can also be 
derived from the finite U result, Eq. ([B]), by taking the 



large U limit {U > e^, i7 > F): 
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(43) 



with A = — Similar procedures as in Appendix [Cl give 
the matrix element i'^'^"^(A) for the dot occupation in 
the infinite U Anderson model in equilibrium to be 



,SBA 



(A) 



2F 



(A - e,)2 + (2F)2 



(44) 



In going to the out-of-equilibrium regime (/ii ^ /Z2) we 
follow the same phenomenological method as for the fi- 
nite U case. The result for the spin-fluctuation and 
charge-fluctuation contributions to the dot occupation 
are given by 



""^^"^^ r V(e,-'A)2+4F2 



2F 



(A - e,)2 + (2F)2 



(45) 



We shall again check the consistency with the exact result 
for the dot occupation in equilibrium, namely 



4/ dA <7,(A);.r^(A) 

ID 



B 



= 4 / dA a,(A)(^^(A)/r (A) + ^.^(A)/,r(A)) . 
Jd 

Here D is related to the bandwidth and B is determined 
by the equilibrium Fermi energy /ii = ^2 = 0. /^(A) 
and /^(A) are expressed as 



/r(A) 



/r(A) 



(A - Bf + (r-)2 
2r 



[A-B- edf + (2r)2 ■ 
Here the Kondo scale used in (A) takes the form^ 
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The results for the dot occupation and Friedel sum rule 
check in the infinite U case are shown in Fig[TTJ Again we 
see a nice match between our phenomenological approach 
and the exact result for Iy'I 7^ ^-^^d some mismatch in 
the mixed valence region] -^j ~ 0. This is consistent with 
the results for finite U. 





FIG. 11: Top: (hd) vs ^ for exact TBA result and pSBA. 
Bottom: Linear response conductance dl /dV\v^o vs y' f^'' 
exact result (FSR) and pSBA in the infinite U Anderson 
model. Y = —100. Similar to the case of finite U the com- 
parison nearby mixed valence region (e^ ~ 0) is poorer. 

The corresponding spin and charge fluctuation matrix 
element for current, J^{A) and J^(A), are given by 



^^(A) = l- 



Cd- A 



V(ed-A)2+4F2 
2F2 



(A-e,)2 + (2F)2 
The current expectation value is given by 

(/) = I dAcT(A)(J^(A)/r(A) + J^(A)/-(A)) 



(46) 




FIG. 12: Js(A) and A (A) vs Bethe momenta A (scaled by F) 
in infinite U Anderson model. ^ = —4 in this graph. Similar 
graph appears for finite U case with x-axis replaced by real 
part of Bethe momenta x{\). 



where Bi and B2 are related to /ii and /i2 by minimizing 
charge free energy F 



F 




dAcr(A)(A-/ii) 



dA(T(A)(A-M2) 



Before we proceed to discuss the numerical results for 
current vs voltage in this infinite U model let us look at 
the structure of (A) and (A) as a function of A as 
shown in Fig. [121 A here represents the bare energy of 
the quasi-particle and plays the same role as x{X) in the 
finite U Anderson model. (A) alone would reproduce 
the main feature in the Friedel sum rule for <C 0. In 
this region the linear response conductance comes mainly 
from the spin fluctuations. The upper plot of Fig. [Infixes 
ed and shows J^(A) vs A. We may also fix A = (in the 
sense of choosing the equilibrium Fermi surface energy at 
A = 0) and plot ^^(e^) vs e^. In this way we can see that 
Jooi^d) VS Erf reproduces the overall structure of the linear 
response conductance from the Kondo region (e^ < 0) to 
the mixed valence regime (e^ ~ 0) . Therefore we identify 



the phase shift 



contributing to J^(A), as the 



phase shift related to spin-fluctuation. 

J;^(A) gives a Lorentz shape in bare energy scale A. 
This structure is akin to the charge fluctuation side peak 
with peak position at energy scale around ej, as seen from 
lower plot of Fig. [T^l Thus we identify the phase shift 
Sp+ + 5p- , contributing to J^(A), as the phase shift re- 
lated to charge-fluctuation. These structures also apply 
to the case of the finite U Anderson model. 

Now let us discuss the out of equilibrium numerical 
results. The voltage is again driven asymmetrically by 
fixing /ii ~ and lowering fi2- The exact dot occupation 
vs voltage for different for infinite U and U = 0, y ~ 
—6 case (black dots) are shown in Fig. [T^l We see again 
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FIG. 13: (fid) vs Y in infinite U Anderson model (for Red, 
Blue, and Purple dots corresponding to ^ = —6, —5, —4. The 
Black dots are (7 = and ^ = — 6 case shown for compari- 
son). Y ~ —100 in this graph. 



the dot occupation decreases slowly at low voltage and 
develops an abrupt drop at a voltage scale corresponding 
to impurity level ed- Also notice the apparent difference 
between the U = plot (black dots) and the C/ — > oo case 
(red dots) and for the same value of For U oo, the 
dot occupation at large voltage is around 0.65 for ^ ^ 
which is consistent with the result of the finite U case 
when Y is large (cf. Section II D). In contrast the non- 
interacting case {U = 0) shows that (n^) — > 0.5 at large 
bias. 

The phenomenological current vs voltage and the cor- 
responding differential conductance vs voltage are plot- 
ted in the top figure of Fig. [HI Again we see the zero 
bias anomaly and a broad charge fluctuation side peak 
in the differential conductance vs voltage. The scaling 
relation of differential conductance vs voltage expected 
in small voltage region can also be extracted by rescaling 
the voltage by T^* as shown in bottom figure of Fig. [TH 
Here T^* is given by 

TT 

Notice this Tj^* differs from with a factor of two 
within the exponent. This factor of two difference rep- 
resents the difference in the curvature of the parabola as 
function of ed (the logarithm of half width at half maxi- 
mum of the Kondo peak vs ed shows parabolic curve as 
in inset of Fig. [5] for finite U case). This factor of two ra- 
tio bears even closer resemblance to the results shown in 
Ref.d. Note that in bottom figure of Fig.[T4]the positions 
of the side peak are different and show no universality in 
that region. It shows universality for < 1. 
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FIG. 14: Top: vs in infinite U Anderson model. Inset 
shows the I—V curves for these parameters, y = —100 in this 
graph. Bottom: ^ vs ^ shows the scaling relation nearby 
zero voltage for ^ = —6, —5, —4 (Blue, Purple, Brown). 



IV. CONCLUDING REMARKS 

In this article we have explicitly computed the non- 
equilibrium transport properties in the Anderson model 
for all voltages using the Scattering Bethe Ansatz. In the 
case of equilibrium we have also shown the equivalence of 
traditional Bethe Ansatz and Scattering Bethe Ansatz by 
evaluating dot occupation in equilibrium. For the expres- 
sion of current we have introduced phenomenological dis- 
tribution functions to set the weight for spin-fluctuation 
and charge- fluctuation contributions to the current. The 
result shows correct scaling relation in Kondo regime as 
well as satisfying the Friedel sum rule for linear response 
for large jr- 

Other interesting quantities, such as the nonequilib- 
rium charge susceptibility or the usual charge suscep- 
tibility, are computed numerically via exact expression 
for dot occupation as a function of voltage and impurity 
level. We believe this is the first report of an exact com- 
putation of the dot occupation out-of-equilibrium and it 
may have interesting application in quantum computing 
as we understand more the dephasing mechanism. We 
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have also compared our results with perturbation calcu- 
lation and experimental measurement of nonlinear differ- 
ential conductance of a quantum dot. 

The major difficulty we encounter by using SBA comes 
from the single particle phase shift for complex momenta 
which leads to a breakdown of steady state condition 
when out of equilibrium. One possible issue resulting 
in this is the local discontinuity at odd channel Sop, 
the choice we made to enable us to construct a scat- 
tering state with fixed particles from lead 1 and lead 
2. It can be proved that without this choice we can- 
not write down fixed number of particles incoming from 
each lead^ in this Anderson impurity model and simi- 
larly for IRLM. The other issue in the study for Ander- 
son model is whether we shall include all possible bound 
states in the ground state construction. From the math- 
ematical structure we shall choose 4 type of bound states 
but the results from charge susceptibility seems to sug- 
gest 2 type of bound states is the correct choice. To check 
whether this is in general correct we plan to come back 
to study the whole spectrum, which include bound state 
when Bethe energy higher than impurity level, of IRLM 
as this model bares structure similarity to the Anderson 
model described in this article. Following the SBA on 
IRLM^ there are lots of numerical approach and different 
exact methods^S. developed for this model and detailed 
comparison for different approaches is desired for better 
understanding its physics and scaling relation. By learn- 
ing how to deal with complex momenta in this model we 
may also find the rule which may lead us to the exact 
expression for current in this Anderson impurity model. 
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Appendix A: Discussion of 2 strings vs 4 strings 

As we have discussed in the main text the bounded 
pair, formed by p^(A) = a;(A) =F iy(A), can be formed by 
quasi- momenta from lead 1 or lead 2. We have shown the 
results for two type of strings (bound states) . Namely the 
strings are formed by {ij} = {11,22} with i, j denoting 
incoming lead indices. In this section we discuss the case 
of 4 type of strings and show thier corresponding numer- 



ical results in out of equilibrium regime (In equilibrium 
the 2 strings and 4 strings give the same result for dot 
occupation). 

The density distribution for the Bethe momenta (ra- 
pidities) is denoted by (Jij{X) with {ij} = {11, 12, 21, 22} 
indicating the incoming electrons from lead i and lead j. 
The crij{X) is given by 

4a.,(A) = -i^- E rdA'if(A-A>.,(A') 

(Al) 

The factor of 4 indicates 4 type of possible configurations 
and the constraint of exclusions in rapidities A in solv- 
ing the quantum inverse scattering problem. The idea is 
that in equilibrium four type of distributions are equally 
possible for each bound state bare energy 2a::(A). The Bij 
play the role of chemical potentials for the Bethe- Ansatz 
momenta and are determined from the physical chemical 
potentials of the two leads, fii, by minimizing the charge 
free energy, 

/•oo 

F = E(^' - ^'^^) = E / (^(^) ~ '^'')'^(^) (^)'^^ 

with CT(i) = 2(711 + o'i2 + (721 the lead 1 particle density 
and (7(2) = 2(722 + o'i2 + o'2i the lead 2 particle density. 
In the case of fii > /i2 we have Bn < B12 = B21 < -B22 
for this finite U Anderson model but the equation for 
(7y(A) is the same for different combination of i and j. 
The reason is we put a quasi-hole state, rather than a 
quasi-particle, in the integral equation Eq. (|Al|) similar 
to the treatment of Wiener-Hopf approach. For example, 
for i?ii < A < ^22 there could be three type of quasi- 
particle state {ij} — {11, 12, 21} and we put {ij} ~ {22} 
state as quasi-holc state. This hole state still count one 
weight of the probability of 4 distributions and therefore 
the factor of 4 on the left hand side of Eg. dAip retains 
even out of equilibrium. Similar idea is also applied in 
two type of bound state (strings) solution. 

Other than their differences in the density distribution 
the computations for the current and dot occupation ex- 
pectation value are quite similar to the two strings case. 
We show their numerical results in the following. 

The differential conductance vs voltage as shown in 
Fig llSl obtained by taking numerical derivative on cur- 
rent vs voltage data, essentially gives the same picture as 
in two strings case, namely a sharp Kondo peak nearby 
V = and a broad side peak corresponding to charge 
fluctuations. In the case of {rid) vs V, however, there is 
an additional feature occurring at an energy scale higher 
than the energy scale of the charge fluctuation side peak 
(corresponding to the voltage position of 2nd peak shown 
in the inset) as shown in Fig. 1161 This is especially ap- 
parent if we looked at the nonequilibrium charge suscep- 
tibility as shown in inset of FigfTEl 

As we do not expect there should be any further charge 
fluctuations, we rule out, by physical argument, the pos- 
sibility of 4 strings configuration. 
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FIG. 15: vs jfr for f/ = 8, r = 0.25 and various ed from 
Ed = — Y to Ed = 1. The inset is the enlarged region nearby 
zero voltage. 
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FIG. 16: (rid) vs ^ for different f/, F = 0.25 and ed = -f. 

The inset is — ^^^le^ "^^ ^ voltage. A third peak shows up 
in U = 4: case. 



Appendix B: Two particles solution and choice of Sop 

For the two particles solution we follow similar con- 
struction in P B Wiegmann and A M Tsvelick's work^ 
and the Scattering Bethe Ansatz approach developed by 
P. Mehta and N. Andrei^. Since Eq.(IT]) is rotational in- 
variant the spin quantum number is conserved. We show 
the solution with both particles with spin singlet incom- 
ing from lead 1 as an example in the following. Spin 
quantum number in z direction Sz is a good quantum 
number and we can write the two particle solution of 



Sz = state as: 

I*) = I J dxidx2{Ag{xi,X2yjl^{xi)il;l^{x2) 
+ Ch{xi,X2)ipl^{xi)ipl^{x2) + Bj{xi,X2){ipl^{xi)ijl^{x2) 
-4M)^Ux2))} + J dx{Ae{x){^l^{x)dl-4^{x)d\) 
+ Bo{x){ijl^{x)dl - i;l^{x)d\)) + Amd\dlj\0) 

Here A, i?, C are arbitrary constants to be determined 
later. To satisfy H\^) = E\'^) = (fc + we have: 

= [-i{dxi + dx^) - E]g{xi,X2) 

+t[S{xi)eix2) + 6{x2)e{xi)] (Bl) 

= [~i{d^,+d^,)-E]h{xuX2) (B2) 

= [-1(8^, + d^,) - E]j{xi,X2) +tS{xi)o{x2)iB3) 

= {-id^- E + ed)e{x)+tg{0,x) + t6{x)m (B4) 

= {-id,- E + ed)o{x)+tj{0,x) (B5) 

= ([/ + 2ed)m + 2te{0) - Em (B6) 

For U — the model becomes non-interacting and the 
two particles solution becomes direct product of two one 
particle solutions. 

= J dxidx2{igkixi)iljl^{xi)+hkixi)'i(jl^{xi)+ekd\S{xi)) 

{gp{x2)ipl^{x2) + hp(x2)tpl^{x2) + epdld{x2)) 
+ (gp{xi)tpl^{xi) + hp(xi)^l^{xi) + epd^^6(xi)) 
{9k(x2)tpl^{x2) + hk{x2)tpl^{x2) + efcrf]^(5(a;2))}|0) 

Therefore at [/ = we have: 



9{xi 


,X2) 


= 9k{xi)gp{x2) + gk[x2)gp{xi) 


h{xi 


,X2) 


= hk{xi)hp{x2) + hk{x2)hp{xi) 


j{xi 


,X2) 


= 9k{xi)hp{x2) + hk{x2)gp{xi) 




e{x) 


= ek9p{x) + epgkix) 




o{x) 


= ekhp{x) + Cphkix) 




m 


— 'ZcpCk 



Now for [/ 7^ we shall derive the solution of this form 

g{xi,X2) = Zkpixi - X2)gkixi)gpix2) 

+ Zkp{x2 - xi)gk{x2)gp{xi) (B7) 

Plug Eq. jBTl into Eq. ((BT1) we get 

e{x) = Zkp{-x)gp{x)ek + Zkp{x)gk{x)ep (B8) 
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Plugging above two results into Eq. (jB4p into Eq. (jB6l) we 

get for m = 2Zf^p{0)ekep we have: 

{-idxZkp{-x))gp{x)ek + {-idxZkp{x))gk{x)ep 
~tZkp{-x)ep5{x)ek 

-tZkp{x)ek5{x)ep + 2tZkp{0)ekep = (B9) 



shown above Paia2 = ^1 thus we have: 



2Zkp{0)ekep 

^ 2t{ZkpiO)gp{0)ek + Zfcp(0)gfc(0)ep) 
p + k-U -2ed 



(BIO) 



and 



Now take Zkp{x) — e 
we get tan((/)fcp) = ^k~p){p+k-u-2e,) 
^fep(O) = klt-u'^X. ^kpiO). Define T = f and 
B{k) = k{k — 2ed — U) as in Ref. [13 we can rewrite 
tan((/)fcp) ^ (B{k}-B{p)) ■ 



From Eq. (jB2|) we can write h{xi,X2) as: 

h{xi,X2) = Zl°p{xi - X2)hk{xi)hp{x2) 

+ Z°k;{x2 ~ xi)hk{x2)hp{xi) (BU) 

with arbitrary Z^p{xi — X2). Now write j{xi,X2) as: 

j{xi,X2) = Zlp(xi - X2)gk{xi)hp{x2) 

+ Zll{x2 - xi)hk{x2)gp{xi) (B12) 



again with Z|°(a;i — X2) undetermined. Plug Eq. (jB12[ ) 
into Eq. ()B3|) we get o{x) is written as: 

o{x) = Z%°p{-x)hp{x)ek + Z%°p{x)hk{x)ep (B13) 

Now if we choose Z^'^(xi — X2) = Z^pixi — X2) and plug 
Eq.dHill) and Eq. ((BT3)) into Eq.dH]) we get: 

(-fc + ed)Zkp{-x)hp{x)ek + {-p + ed)Zkp{x)hk{x)ep 
+t{Zkp{-x)hp{x)gk{0) + Zkp{x)hk{x)gp{0)) 
+{-i){dxZkp{-x))hp{x)ek + (-i)(dxZkp(x))hk{x)ep 
= -2 sin(0fep)(/ip(O)efc - hk{0)ep) = (B14) 



To satisfy Eq. (jB14l) we can set hp{Q) — for arbitrary p. 
This can be done by choosing Sop — —4 in Eq.®. Now 
since Z'^°{xi — X2) is arbitrary we can choose Z'^°{xi — 
X2) = Zkp{xi — X2)- Also from Ea. (|B10[) we have 



ZkpiO) 



p + k-2ed 
p + k-U -2ed 



Zkp{0) 



(B15) 



Since the Hamiltonian in Eq.([T]) has rotational invariance 
the general form of scattering matrix for particles with 
momentum k,p and spins 01,02 is given by: 



Salal{k,p) = b{k,p) + c{k,p)Palai 



(B16) 



where Pala^ = 5 (1^1 • laa +'?aj ' (Taa ) IS the permutation 
operator in spins. For antiparallel spins (singlet state) as 



b{k,p) - c{k,p) 



Zkpjx > 0) 
Zkp{x < 0) 
B{k) - B{p) - i2UV 
B{k) - B{p) + i2UV 



(B17) 



For the triplet state {Pala2 = 1) the interaction term 
with the impurity is absent and the particles passing 
through each other without changing their phase 



b{k,p) + c{k,p) = 1 



(B18) 



Thus from Eq. (|B17p and Eq. (jBlSp we get the two particle 
S-matrix as: 

sik vt^"'^ - (i3(fc)-i?w)i:::j+»2[/rps:i 

^^'''P^"""^ - B{k)-B{p)+i2UT ^^^^^ 

Thus the integrability of two lead with Anderson type 
dot system is the similar to the integrability of one lead 
Anderson model. 

The choice of identical two particles S-matrices (by 
choosing Sop = —4) enables us to construct the scat- 
tering state labeled by lead indices by choosing appro- 
priate B, C in this even-odd basis. For example, if 
both particles are coming from lead 1, we shall choose 
{A,B,C) = ^o(tt, r?-, 72-) such that the amplitude of 

^2 ''1''2 ^1 

incoming state from lead 2 is zero {Aq being an overall 
renormalization constant). We can therefore label the 
eigenstate by the incoming state from lead i and/or lead 
j. Without this Sop term we cannot write back from even- 
odd basis to lead indices basis in this two leads Anderson 
model and similarly in IRLM in Ref. i7|. 



Appendix C: Equivalence of TBA and SB A in 
equilibrium 

Eq. p^ can be proved to be exact by comparing 
with the traditional Bethe Ansatz where d\da) — 
2 dA(Tiinp(A) with impurity density crimp (A) given by 

^imp(A) ^ - / dA'if(A-A')aimp(A') (CI) 

Jb 

The driving term (first term) of Eq. (ICip is expressed by 
bare phase shift dp+ +6p- and thus we can view (Tiinp(A) as 
the dressed phase shift across the impurity. By compar- 
ing Eq. (|Cl|) and Eq.® in equilibrium (tTi(A) = ab{X) de- 
scribing bulk quasi-particle density when Bi ^ B2 ~ B.) 
we get 



dA.imp(A)(^^) 



TT dA 



(C2) 
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by noting that the integration kernel K{\ — A') is sym- 
metric in A and A'. Since the equahty is true for arbitrary 
B we can also rewrite Eq. (jC2|) as 



dA crimp(A) 



dA 



2 / dXatiXy'^iX) 

IB 



and the resulting is given by 



^(A)f^ + y+(A) 
52(A) + y2 (A) 



(C3) 



Now let us show the computation for v^^'^{X). First we 
write one particle state of Eq.((T]) in even channel (with 
Sek = for the moment) as 

\Ka)= j e^^-al^Jx)dx\Q) (C4) 

Solving H\k,a) — k\k,a) we get 

+ Ak) + Bkt = 
l + Ak 



^dBk + 1- 



kBk 



Thus we get A^ 



-3r and Bk ~ ^ — -tt ■ We may 

and efe = Bk to have 



Aki 



also define gk(x) = e 
easier comparison with Wiegmann and Tsvelick's work's. 
The two particles state is obtained by constructing prod- 
uct of two Qfjp o- {x) particles state with appropriate two 
particles S-matrix expressed in Zj^+f^~ (xi — X2)- 

In principle we shall use , Ni, N2) as the many body 
state to compute expectation value. However the simpli- 
fication here, similar to the case of IRLM in Ref.?, is that 
different A (corresponding to different p{X)) are orthogo- 
nal to each other in L — cx) limit. Thus the many body 
expectation value can be obtained via two body compu- 
tation and the rest just get canceled by normalization 
factor. To put it more explicitly let us denote pi as the 
real part of the complex pair pf . Different \pi) is orthog- 
onal to each other under the condition of size of the leads 



taken to infinity, or 



{pi\pi 



— )■ as L — >■ 00 for i 7^ j. Thus 



the evaluation of matrix element for operator 6 is given 

by 



{pi,P2, . . . |6bl,P2, 
(j>l,P2,-.-\pi,P2,- 



E 



{pi\o\Pi 

ipl\pt) 



. Based on this result we demonstrate the explicit compu- 
tation for dot occupation by two particles wavefunctions 
in the following. 

Denote | VP) as the two particles solution. We may write 
spin singlet state as 



I*) 



dxidx2A[e''^''^'+P^'^Zkj,{x^ - X2)alf^^^{xi)al^^^{x2)] |0) 

dxidx2\^Zkp{xi - X2){gk{xi)gp{x2)il^\(xi)i)l^^{x2) + gk{xi)ep'4)\{xi)d\5{x2) 
+ ek9p{x2)d].5{xi)il)\{x2) + ekepd].d''^5{xi)5{x2)} - Zkp{x2 - Xi){gk{x2)gp{xi)'4)l^{x2)^l^{xi] 
+ gk{x2)epi'li^{x2)d\.5{xi) + ekgp{xi)d^^^5{x2)i'l^{xi) + efcepdjd|j(a;i)(5(x2)}||0) 

dxidx2[Zkp{xi - X2)gk{xi)gp{x2) + Zkp{x2 - Xi)gk{x2)gp{xi)]il;l^{xi)ipl^{x2) 

dx[Zkp{x)gk{x)ep + Zkp{-x)gp{x)ek]{ipl^{x)dl - ipl^{x)d\.) + 2efcepZfep(0)d|(i|||0) 



With A denoting anti-symmetrization and Zkp{0) 

k+p-2ea y ,r.\ 

k+p-U~2ea ^kp\^)- 



Solving H\k,a;p,—a) — {k + p)\k,a;p, —a) we obtain 



Zkpixi - X2) = 0{xi - X2) + 



{k~p){k + p- 2ed -U)- iUe 
{k-p){k + p- 2ed ~U)+ iUt^ 



0{X2 - Xi) 
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For the case of bound state the two particle S-matrix is 
given by Zj^+j.- {xi — X2) = 0{xi — X2) = 0^2- The nor- 

I 



a 



Note that the even channel bound state can be written as 
sum over bound state of {11, 12, 21, 22} (4 strings type) 
or {11,22} (2 strings type) with the same real part of 
energy k = x{X). This can be viewed as the consistency 

I 



j dxidx2 Ol2\gk+{xi)gk-{x2)\'^ 



j dxe{x)\gk+ (a;)efc-p 



j dxe{-x)\gk-{x)ek+\'^ 



with ^fe+fc-(0) = 2(k-.f-u ^k+k-{Q) and Zk+k-{Q) = | 
based on our regularization scheme. By expressing k = 



malization factor and matrix clement of dot occupation 
given by the even channel two particles wavefunction are 



I 

counting from Fock basis to Bethe basis as electrons in 
lead 1 and lead 2 has 4 fold degeneracies in its initial 
state (2 different spins in each lead). Also note that 



I 

x{X) and = J/(A) and taking L ^ 00 thus preserving 
4- terms only we get 



= j dyidy2 j dxidx2{0f29k+ {yi)9k- {y2) + Oligk+{y2)9k- ivi))* 
^{0i29k+{xi)9k-{x2) + 02i9k+{x2)9k- (a;i))^(a;i -y\)5{x2 -2/2) 

+2 j dy j dx[0{y)gk+iy)ek- +0{-y)gk-{y)ek+]*[9{x)gk+{x)ek- +0{-x)gk-{x)ek+]S{x-y) 

+4(efe+efe- Zk+k- (0))* (efe+ e^- Z^+k- (0)) 
^ j j (^xi^iy)9k+iy)ek- + 0{-y)9k-{y)ek+]*[0{x)gk+{3:)ek- + 0{-x)gk- {x)ek+]6{x - y) 

+8(efe+efc- Zk+k- {0))*{ek+ek-Zk+k- (0)) 
= 2! / dx[e{x)\gk+{x)ek-f + e{-x)\gk-{x)ek+f]+Mek+ek-Zk+k-{0)f\ 



J dXidX2\e'^''^ ""'H^l + 01 Ak+){e2 + 02 Ak-)\%2 

dxidX2 e-^^''^'''-''''>\0i02012 + 0102012 Ak+ + 0102012 Ak+Ak-\^ 



L 1 - e-^Cfci 

{k+iik)x(o(_„ 



1 - e-2«fc^ 



2a 



\Ak+\ 



{l + \A,+ A,-f) + 
J rfa;6l(x)|e'('=+'«'=)^(6'(-x)+Afc+6'(x))efe-|2 = ^ dxe-2«'=^|Afc+efc 

k- cd + i^k- «r 



1 
26 



t 



1 
26 



t 



k- ed + i^k+ «r 



k - €d + i^k k - ed - i^k + 

dx0{-x)\e'^''-'^''^''{0{-x) + Ak-0{x))ek+\^ = J dxe2«'=^|Afe-efc+| 



1 

26 



t 



k-€d + i^k + «r 
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16y(A)r2 / i(A) 




i^{\) + yl{\) (i2(A)+e(A))(i2(A)+y2(A)) \2i{\)-U 



By expressing j/-^^'^(A) and u^^^{X) in A explicitly we 
see that i/^^^(A) = i/^^^(A). Since {^2^4^^} = 
2 dAcrjmp(A) in TBA we have proved that the expecta- 



tion value evaluated by the state wc constructed is exact 
and the equivalence of SEA and TBA in equilibrium in 
this two-lead Anderson model. 
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